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■ An exact expression for the gravitational field strength in a self-gravitating dust 
ff^ , continuum is derived within the Lagrangian picture of continuum mechanics. From 

^ \ the Euler-Newton system a transport equation for the gravitational field strength is 
formulated and then integrated along trajectories of continuum elements. The result- 

. ing integral solves one of the Lagrangian equations of the corresponding Lagrange- 

\ Newton system in general. Relations to known exact solutions without symmetry in 

ly-^ I Newtonian gravity are discussed. The presented integral may be employed to access 

^) • the non-perturbative regime of structure formation in Newtonian cosmology, and 

CJ ! to apply iterative Lagrangian schemes to solve the Lagrange-Newton system. 
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1 The problem 



Let us characterize the state of a self-gravitating "dust continuum" (i.e., a 
continuum without pressure) at the initial time tj by a velocity and a density 
field v(x, ti) and ^(x, t,). As usual, the fields are represented in a non-rotating 
Eulerian coordinate system x. We are interested in the evolution of the contin- 
uum governed by the Eulerian evolution equations for the velocity field v(x, t) 
and the density field ^(x, t): 
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The evolution of the continuum is constrained by Newton's field equations 
for the gravitational field strength g(x, t), (which is equal to the acceleration 
field according to Einstein's equivalence principle of inertial and gravitational 
mass) : 



where A is the cosmological constant, and the density field obeys ^ > 0. 
(The introduction of a gravitational potential $ defined by g =: — V$ (the 
existence of which is guaranteed by (4)) will not be needed.) Eqs. (1 - 4) form 
the Euler-Newton system of equations [18], [19], [20], [15]. 

Integration of this system is relevant in the context of Newtonian cosmology, 
where the "dust continuum" models the (together with the source A) dom- 
inant collisionless "dark matter" and "dark energy" in the Universe, and is 
mostly studied perturbatively, i.e. perturbative solutions at a homogeneous- 
isotropic solution of the Euler-Newton system are found in the Eulerian and 
Lagrangian representations of the system (see, e.g., [22], [23], [16], [3] and ref- 
erences therein). Besides these analytical schemes the problem is commonly 
treated by employing a variety of numerical integrators [4]. Exact solutions 
for inhomogcncous fields are known for special spatial symmetries like planar 
and spherical symmetry; three-dimensional solutions without symmetries are 
rare and will be put into perspective in Sect. 5. 

We proceed as follows. Wc start by deriving a transport equation for the grav- 
itational field strength, working in the traditionally more emphasized Eulerian 
picture (Sect. 2). Then, we move to the Lagrangian picture and integrate the 
transport equation along the flow-lines of continuum elements, flrst by setting 
the cosmological constant A equal to zero, in Sect. 3. The result is generahzed 
to the case of a non-vanishing cosmological constant in Sect. 4. Sect. 5 provides 
relations to known exact solutions and presents an iteration scheme to solve 
for the trajectory field. Sect. 6 gives an outlook on applications of the presented 
integrals in Newtonian cosmology. 



V-g = A-47rG'^ ; 
V xg = , 



(3) 
(4) 
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2 Transport equation for the gravitational field strength 



We first recall some equations that have been obtained earlier [9] , and which 
provide an alternative formulation of the Euler-Ncwton system: we can com- 
bine the evolution equation (2) with Eq. (3) to obtain the following evolution 
equation for the gravitational field strength ([9] Eq. (7c*)): 

^g - Av = (v • V)g - v( V • g) + V X T , (5) 

where ^ := + v • V is the total (Lagrangian) time-derivative, abbrevi- 
ated sometimes by an overdot. (Hereafter, we shall drop the subscript |x for 
notational ease.) 

The vector field r specifies the freedom after formally integrating the diver- 
gence of equation (5), and can be interpreted as the vector potential of the 
current density j = ([9] Eq. (15a*)): 

47rG'j = -V^* - V X r , (6) 

where $ is the gravitational potential as introduced above. Eq. (6) is equivalent 
to Eq. (5) using the definition ^g = ^g + v • Vg and the source equation for 
the gravitational field strength (3). 



In passing 1 mention that the field r plays the role of a "magnetic" field 
strength in a Maxwell-type analogy, while g corresponds to the "electric" 
field strength; we may cast Newton's field equations together with Eq. (6) 
into the familiar form: 

V ■ g = A - 4:TiGg V • r = ; 

Vxg = VxT = ^g- 47rG'j . (7) 

The constraint V • r = is a gauge condition (transverse gauge) imposed on 
the vector potential r. It is always possible to require this condition according 
to the following reasoning: the curl of the vector field r remains unchanged, 

if we add to r the gradient of an arbitrary scalar field, r' := r + 'Vip; the 
divergence of r itself is not specified by the Euler-Newton system, so we can 
find ijj such that V ■ r = from the following equation: 

AV'^V-r' . (8) 

There always exist solutions to Eq. (8) according to a theorem by Brelot (see, 
e.g. , [17] ) , if the spatial average of the source of Poisson's equation (8) vanishes. 
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Using this gauge condition, we find from (6) the simple relationship: 

AtiG V X j = - V X (V X t) = At - V (V • t) = At . (9) 



Let us return to the derivation of the transport equation. We employ the 
following vector identity: 

V X (g X v) = (v ■ V)g - v(V ■ g) - [ (g ■ V)v - g(V • v) ] . (10) 

With this identity, Eq. (5) can be cast into a transport equation for the grav- 
itational field strength: computing g-^ (g/f?), using (5), and inserting the con- 
tinuity equation (2) in the form g/g = —V • v, we obtain: 

T- (- I = (- ■ I v + -v + -V X f ; f:=T + gXv, (11) 

dt \gj \g J g g 

which is the key-equation of the present work. Eq. (11) (for A = 0) is formally 
similar to Beltrami's transport equation for the vorticity u) :— | V x v reviewed 
in APPENDIX A. 

In the next section we shall integrate the transport equation (11) after intro- 
ducing the Lagrangian form of the Euler-Newton system. 



3 Lagrangian integral of the transport equation for A= 

The Lagrangian description is based on integral curves x = f(X, i) of the 
velocity field v(x, t): 

^=v(f,t) ; f(X,t,)=:X . (12) 

In other words, f is the position vector field that locates a fluid element, 

indexed by its initial position vector (the Lagrangian coordinates Xj), in Eu- 
lerian space at a given time t. Introducing this family of trajectories, we can 
express all Eulerian fields, e.g., the velocity v, the gravitational field strength 
g, the density g, and the vorticity oj {cf. APPENDIX A) in terms of the field 
of trajectories x = f (X, t) as follows: 



v = f(X,i) ; g = f(X,t) ; (13) 
^=ft(X) J-i ; (14) 
w = • Vof , (15) 
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with the Jacobian of the transformation from Eulerian to Lagrangian coordi- 
nates J := det(/j|fc(X, t)) > 0, oJi := a;(X, tj), and the nabla-operator with 
respect to Lagrangian coordinates denoted by Vq- (Throughout this paper a 
vertical slash | denotes partial differentiation with respect to Lagrangian coor- 
dinates, which commutes with the Lagrangian time-derivative, while Eulerian 
spatial differentiation is abbreviated by a comma.) 

Eqs. (14) - (15) are the known Lagrangian integrals of the Euler-Newton 
system, i.e., they represent an Eulerian field as a functional of f. As a result 
of the definitions (13) and the density integral (14) both Eulerian evolution 
equations are solved exactly for any given trajectory field in the Lagrangian 
picture. Transforming those fields back to Eulerian space wc need that the 
transformation f is invertible, i.e. J > 0, defining regular solutions (for more 
details the reader may consult the review [16]). 

While the Eulerian evolution equations are represented by definitions and 
Lagrangian integrals, the relevant equations in the Lagrangian picture are the 
Eulerian field equations (3) and (4), which are transformed into a system of 
Lagrangian evolution equations by virtue of the following transformation of 
the field strength gradient gij — -^gi- 

9i,j — gi\khk,j — -;^^klm^jpqgi\kfp\ifq\m ] (16) 

the gradient of the inverse transformation from Lagrangian to Eulerian co- 
ordinates, h = (which we require to exist), was expressed in terms of f 
through the algebraic relationship 

^ij ~ '^ij^ ~ 0'd{Jij)J ^ = -^^ike^jmnfmlkfrili ■ (17) 

In view of (3) and (4) we obtain with (16) the following set of four Lagrangian 
equations ([8] (A = 0), and [9] (A 7^ 0); i,j,k — 1,2,3 with cyclic ordering; 
summation over repeated indices is understood): 



Note that these equations only involve the dynamical variable f by virtue 
of g = f . Inserting this definition into the equations above, we obain a set of 
four evolution equations furnishing the Lagrange-Newton system. (Alternative 
forms of these equations may be found in [12] and [16].) 



We now integrate the transport equation (11) along trajectories f of contin- 
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uum elements indexed by X. The system (18,19) is a set of nonlinear partial 
differential equations coupled in a complicated way, and certainly more in- 
volved than their Eulerian counterparts that are linear field equations. (Of 
course, the complete Eulerian system is nonhnear due to the couphng of the 
field strength to the density, which itself has to obey a nonlinear equation.) 
With regard to the complicated nature of the Lagrangian equations it may be 
surprising that one can find - as we shall see - a relatively simple expression 
for g as a functional of f, extending the set of known integrals Eqs. (14) - 
(15). With the condition 

V X f = (20) 

we single out a class of motions that admits a quasi-local Lagrangian integral. 
The term 'quasi-local' is to be understood as follows: if we expect g = f to 
be a functional of f, then the trajectory field is represented locally, while the 
initial data are still constructed non-locally according to the structure of the 
theory, i.e. we have to prepare initial data such that all the field equations 
(3,4) are satisfied initially. However, we also expect a general integral, if it 
exists, to contain non-local parts that arise when the imposed restriction (20) 
is relaxed. Requiring (20) to hold is sufficient to enable us to integrate Eq. (11), 
but it is not necessary for the following Proposition 1. 

Setting first A = we make the following ansatz: 

^ = (s-Vo)f ; s = s(X,t) . (21) 
Q 

Applying the total time derivative to this ansatz, and using the identity 

(s-Vo)f =((s-Vo)f)-Vv , (22) 
we compare the resulting terms with the transport equation (11) to obtain: 

Q 

(s • Vo)f = , i.e. , in general : s = const. = — ; G := g(x, ij) . (23) 

Qi 

Using the integral for the density field, Eq. (14), we arrive at an exact integral 
for the gravitational field strength (this result has been first given in [11], 
however, without providing the derivation): 



Proposition 1 (Lagrangian Integral for A = 0) 
The Lagrangian integral 

g,^(G^Vo)f . j^det(/,|,) (24) 

solves the Lagrangian equation (18) for A = exactly, while the remaining 
equations (19) are reduced to constraint equations restricting f. 
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(The proof follows by explicitly inserting the integral into the Lagrangian 
equations, either manually or using an algebraic manipulation system, e.g. 
REDUCE. We outhne an analytical proof in APPENDIX B.) 

Note that this integral has its counterpart in the classical integral (15) for the 
vorticity as reviewed in APPENDIX A. 



4 The case of a non— vanishing cosmological constant 



The introduction of a cosmological constant in the Euler-Newton system cor- 
responds to the introduction of a background field. In Newtonian cosmology 
such a background is commonly installed in terms of a "Hubble flow", i.e., a 
homogeneous-isotropic solution of the Euler-Newton system that is furnished 
by the deformation field: 

fn = a(t)X ^ = -X ; = QH(ti)a~^ , (25) 

a 

with the "scale-factor" a{t), the linear "Hubble velocity" vh, and the ho- 
mogeneous background density Qnit)- Inserting and gn '■= into the 
Lagrange-Newton system (18), (19), we obtain the well-known cosmological 
equation: 

3- ^ A-AttGqh ■ (26) 
a 

Together with as given in terms of a{t) in (25), this equation determines the 
solutions a{t). Deviation fields from the "Hubble fiow" are then introduced. 
(Such a description hes at the basis of applications to cosmological structure 
formation models, cf. Sect. 6.) In this line, the cosmological constant plays the 
role of a static background, which introduces qualitatively new features into 
the solutions. 

Relaxing the restriction to a vanishing cosmological constant, we have to gen- 
eralize the ansatz (21): since A is constant in space and time, we expect a 
deviation from the field strength in (21) that is proportional to the displace- 
ment defined by f: 

K- -Df 

= (s • Vo)f , (27) 

Q 

with D — const. Applying the total time derivative to the new ansatz (27), 
using again the identity (22) and comparing the resulting terms with the 
transport equation (11) we first obtain: 



D 
Q 



f(V-v) - (f • V)v + - 

Q 



A- . (28) 
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With the following identity we are able to evaluate the constant D: 

V • [ f (V • v) - (f • V)v ] = V • [ v(V ■ f ) - (v ■ V)f ] = V • (d - l)v , (29) 

where the dimension d of the continuum has to be explicitly taken into account. 

Since all equations are thus far understood up to a vector potential that we 
neglected, we obtain in view of (28) and (29): D = A/d. Using the Lagrangian 
integral for the density (14), we finally arrive at a more general exact integral 
for the gravitational field strength: 



Proposition 2 (Lagrangian Integral for A 7^ 0) 

In d dimensions of the continuum the Lagrangian integral 

g,,^ (C_Vo)f ^A^ ; C:=G-^X ; J = det(/.|,) (30) 

solves the Lagrangian equation (18) exactly, while the remaining equations 
(19) are reduced to constraint equations restricting f. 

(Again, the proof follows by explicitly inserting the integral into the La- 
grangian equations, either manually or using an algebraic manipulation sys- 
tem, e.g. REDUCE, see APPENDIX B for an analytical proof.) 

The exphcit dependence of the integral (30) on the dimension d of the contin- 
uum can be easily understood by noting that from ^ = 0, V • g = A, we have 
the integral g = which is confirmed by calculating the divergence of this 
integral with respect to x = f (X, t). 



5 Relation to exact solutions and definition of an iteration scheme 



Let us first consider the integral (24) for A = 0. It provides the (with respect 
to (18) without any restriction) exact field strength along any given family 
of trajectories. Inserting the integral into the remaining Lagrangian equations 
(19) delivers constraints that have to be imposed on the given family of trajec- 
tories yielding in particular restrictions on the spatial coefficient functions. (It 
is possible to find a corresponding exact integral for the Lagrangian equations 
corresponding to V x g = which, again, requires constraints by inserting 
it into (18); this and also a formally general integral will be given elsewhere 
within a more general context; the integral of (18) is favoured with regard to 
applications, cf. Sect. 6.) 
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Alternatively, the integral (24) may be regarded as a set of three partial dif- 
ferential equations for the components of the trajectory field by virtue of the 
definition of g (13). Integrating the integral twice with respect to the time 
and fixing the two integration constants by the initial data: f(x, tj) = X, 
v(x, ti) = V(X), we obtain an integral equation for the trajectory field: 

f = X + y{X){t - ti) + f dt' r dt' (G(X) • Vo)f (X, t')J{X, t'y . (31) 

Jti Jti 

This equation suggests, as a manifest structural property of the integral (24), 
the following iterative definition of the trajectory field f: 

f(n+i) ^x+v(x)(t-ti)+ f'dt' f dt' (G(x)- Vo)f("Hx,0'^^"Hx,0"' , 

Jti Jti 

(32) 

with J("):=det(/(^)(X,0-^ , 

where f is the n^^ iterate of the integral. A trajectory field that is inserted 
on the r.-h.-s. of (32) does not in general yield the same trajectory field f on 
the l.-h.-s. of (32). This way we only consider f as a solution of the Lagrangian 
equation (18), iff f is, in the language of dynamical systems theory, a fix-point 
of this iteration, i.e., 

f (n+l) ^ f (n) ^ ^33^ 

for all n larger than some n*. 

Let us now consider some known cases of exact solutions in Newtonian grav- 
ity. Among the simplest of exact solutions we have already introduced the 
homogeneous-isotropic deformation fields (25) and (26). Inserting this so- 
lution for A = into the r.-h.-s. of (32), we recover the same solution on the 
l.-h.-s. of (32) (with careful handling of the integration constants). Thus, in 
this case, the iterated solution is identical to the solution fed into the integral. 

We face the same situation in the case of the exact general one-dimensional 
solution [21], [8]: 

/i = + V,{X,){t - ti) + Gi(Xi)^^^^ . (34) 

Also, the three-dimensional generalization of this solution [8], admitting lo- 
cally (at each X different) one-dimensional motions and imposing no global 
symmetry restrictions, forms a fix-point of this iteration: 

f = X + V(X)(i-i,) + G(X)^^^^ . (35) 

Here, we have to make explicit use of the constraints given in [8] and the iden- 
tity (10) to verify this. (The solution class is defined by three-dimensional 
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initial data that are composed of 2-surfaces with vanishing Gaussian curva- 
ture; the plane-symmetric case is contained in a subclass formed by cylindrical 
2-surfaces.) The basic assumption that was used to derive this class of exact so- 
lutions was the constancy of the gravitational field strength along the solution 
curves: g(X, = G(X). The integral (24) can be viewed as a generalization 
of this assumption, where g changes along solution curves according to their 
directional derivative with respect to the initial field G, weighted in addition 
by the local density. It is interesting that we obtain the solution curves (35) 
as the first iterate of the integral (24), if we start with the trivial trajectories 
f(o) — X. Any further iterate is identical to f(^\ if the constraints quoted 
above are respected. 

Of course, we may start this iteration scheme with generic initial data in order 
to approximate solutions of the Lagrangian system. Iteration will produce 
another approximation that will contain, at each iteration step, higher spatial 

derivatives of the initial data. At this stage it is premature to expect that 
such an iteration scheme may converge in the sense of producing more refined 
models which approach a solution of the Lagrange-Newton system. 



Relaxing the restriction A = 0, we accordingly obtain the more general tra- 
jectory field, which we may write as follows: 



x+y{x){t-ti)+ f dt' f dt' 

Jti Jti 



(G(X)-^X)-Vof J-^ + ^f 

a a 



(36) 



Viewing the above equation as an iteration scheme, then one possible choice 
is to write: 

f(n-M) Af(nH-l) (G(X) - A X) • Vpf ^ 

If we start with a homogeneous-isotropic deformation, f*^°^ = a(t)X, then we 
find a fix-point corresponding to all solutions of Eqs. (25), (26) for d — 2). 

On the other hand, if we start with the trivial trajectories f = X, then the 
first iterate obeys: 

f(i)_^f{i)^G(X)-^X , (38) 

(Jj CL 

admitting the general integral-curves: 



f«=X + A(G(X) + ^^V(X) 



eVi - 1 
e-v^(*-*0 - 1 



(39) 



For the case d — 3 we propose a more refined choice of the iteration scheme 
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(37), motivated by the idea that the first iterate then provides an exact solution 
to the Lagrange-Newton system (18), (19). We write: 

f (n+l) _ Af(-+1) - (^(^) - t ^) • _ ^f(n) 

J(fH) 3 ■ ^ ^ 

With this choice the Friedmannian solutions still form a fix-point correspond- 
ing to the iteration of a homogeneous-isotropic deformation (as above). The 
advantage of this latter choice becomes obvious by iteration of the trivial 
trajectories f = X, since the first iterate then obeys: 

f(i) _ AfW = G(X) - AX , (41) 

admitting the general integral-curves: 



f = X + ^ (G(X) + VA V(X)) [e^ - 

+ ^ (G(X) - VA V(X)) [e-^(*-*^) - l] . (42) 

The trajectory field (42) provides a class of three-dimensional exact solutions 
without global symmetry restrictions that is contained in a subclass of solu- 
tions investigated in [9], which will not be further discussed here (see also [5]; 
for solutions including a cosmological constant see [6], [14], and for related 
rotational solutions [25]). Explicitly, this solution class has been discussed by 
Barrow & Gotz [2] in the context of "no-hair" theorems. 

Although the basic equations tend continuously, in the limit A ^ 0, to their 
restricted set for A = 0, the solution (35) given above can only be obtained 
from (42) as an approximate limit by expanding the exponentials. This shows 
that the presence of a background entails a qualitative difference. 



6 Applications in Newtonian cosmology 



In order to appreciate the application side of the presented integrals, let us give 
a concrete example in the context of structure formation models in Newtonian 
cosmology, which requires the presence of a background. 

We employ the integral (30) for d — 3. Since this integral is exact, we are 
entitled to apply the transformation to a "comoving trajectory field", q 
F(X, t) = :K/a{t), and to introduce deviation-fields from the background so- 
lution: since we can adopt the same Lagrangian coordinates Xj in both cases 
in view of x = a(t)q, a{ti) = 1, F(X, ti) = X, and since the background 
field strength is evolving as gH = oq, we make the ansatz g = aF -|- w(X, t), 
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G = a(tj)X4-W(X), use Friedmann's equation (26), and obtain from (30) the 
following transformed integral for the "peculiar-gravitational field strength" 
w in terms of F : 



C = W — ^n^Mlhly^ ig the integration constant. The iteration scheme is 
defined, for d = 3, according to the choice (40): 

-w["+^l - 47rG^^Ft'^+i] = - ^!^fN , (44) 

with iw["+^l = Y^ri+\\ _^ 2ifF["+^l. The idea is to start this iteration scheme 

a 

for generic initial data defining an approximation to the Lagrange-Newton 
system. Starting the iteration (44) with the trivial trajectory field F^^^ = X, 
which corresponds to a "Hubble flow" , we obtain the equation governing the 
longitudinal first-order Lagrangian perturbation solution in Newtonian cos- 
mology [10], Ft-*^! = X + Pt^l, with p[^l given for all Priedmannian backgrounds 
in [6]; a subclass of this trajectory field is the celebrated "Zel'dovich approxi- 
mation" [26], which is already very successful for the description of structure 
formation in cosmology in comparison with N-body simulations (see the ref- 
erences in the review papers [7], [23], [16]). Any further iterate produces a 
non-perturbative approximation. A numerical implementation of this itera- 
tion scheme is currently under investigation. 



A flnal word on the range of validity of the presented integrals is in order. 
Although we discussed all equations and results in the framework of a "dust 
continuum". Equations (11), (18), (19) and their integrals for g are valid 
for a much wider range of continua, e.g. including pressure forces. In those 
more general cases the representation of g in terms of the trajectory field 
changes, and is no longer given by the simple relationship g = f as in the 
"dust continuum" (see, e.g., [1] for the case of an isotropic pressure). This 
more general context and all details relevant to applications in Newtonian 
cosmology will be presented in a forthcoming paper [13] . 
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APPENDIX A: Lagrangian integral for the vorticity 



We here review a classical result that displays formal similarities to the inte- 
gration procedures used for the derivation of the transport equation for the 
gravitational field strength. The transport equation (11) is formally similar to 
Beltrami's transport equation for the vorticity a; := x v (see, e.g., [24] 
and [10]): 

For curl-free forces, which is the case for Newtonian flows, this reduces to: 

(46) 

There is also a classical integral of Beltrami's transport equation due to 
Cauchy in the case of curl-free forces (see [24] and [10]): 

^^uvVof . ^.^^^x,i,) . (47) 





APPENDIX B: Proof of Propositions 1 and 2 



A simple way to prove Proposition 2 and, being a subcase. Proposition 1 is to 
write the divergence of g as gi^i — gi\khk,i (first without cxpficitly writing out 
h). Computing g^^ from the integral (30), gi — Csfi\sJ~^ + ^ fi, we first have: 



gi\k = 



J 



Cs\kfi\s + Csfi\sk J~Csfi\s 



+ ^fMk 



Multiplying by hk,i and using fi\shk,i = Sgk, we obtain: 



9i,i 



J 



Ck\k + Csfi\skhk,i 'j'^^ 



+ (>kk 



(48) 



(49) 



k — ^ Xk, G 



k\k 



We have to prove that gi^i = A — AnGQiJ~^ . Since Ck = G 
A — AnGQi, X/jifc = 5kk = the first term in Eq. (49) together with 5kk = A, 
already provides the whole equation, gfj j = A — AttGq. It remains to prove 
that, for J 7^ 0, the other two terms in the brackets of Eq. (49) cancel out: 



Gsfi\skhk,i — ■ 



(50) 

fp\rfq\o, and also 

the Jacobian, J = ^eabc^def fa\dfb\efc\f to obtain the equivalent requirement: 



We now write the inverse matrix explicitly, hs 



G s^ipq^sro 



fi\skfp\rfq\o G^ 6ai,Q6^gf 

{fa\dfb\efc\f).j^ ■ (51) 
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We appreciate that, by evaluating the r.-h.-s., we get three identical expres- 
sions, so that by relabelling the summation indices we conclude that the above 
equation holds. ■ 

The REDUCE codes to compute the Lagrangian evolution equations and the 
presented integrals can be obtained from the author. 

References 

[I] S. Adler, T. Buchert, Astron. Astrophys. 343 (1999) 317. 

[2] J.D. Barrow, G. Gotz, Class. Quant. Grav. 6 (1989) 1253. 

[3] F. Bernardeau, S. Colombi, E. Gaztanaga, R. Scoccimarro, Phys. Rep. 367 
(2002) 1. 

[4] E. Bertschinger, Ann. Rev. of Astron. and Astrophys. 36 (1998) 599. 

[5] S. Bildhauer, T. Buchert, Prog. Theor. Phys. 86 (1991) 653. 

[6] S. Bildhauer, T. Buchert, M. Kasai, Astron. Astrophys. 263 (1992) 23. 

[7] F.R. Bouchet, S. Colombi, E. Hivon, R. Juszkiewicz, Astron. Astrophys. 296 
(1995) 575. 

[8] T. Buchert T., G. Gotz, J. Math. Phys. 28 (1987) 2714. 

[9] T. Buchert, Astron. Astrophys. 223 (1989) 9. 

[10] T. Buchert, Mon. Not. Roy Astron. Soc. 254 (1992) 729. 

[II] T. Buchert T., Inhomogeneous Newtonian Cosmogony, Habilitation Thesis, 
Ludwig-Maximilians-Universitat Miinchen (1993); link on ADS. 

[12] T. Buchert, in: Proc. International School Enrico Fermi, Course CXXXII 
(Dark Matter in the Universe), Varenna (Italy), S. Bonometto, J.R. Primack, 
A. Provenzale (eds.), lOS Press Amsterdam (1996) 543. 

[13] T. Buchert, Astron. Astrophys., submitted.; |astro-ph/0601513 

[14] A.D. Chernin, D.I. Nagirner, S.V. Starikova, Astron. Astrophys. 399 (2003) 
19. 

[15] J. Ehlers, in: Grundlagenprobleme der modernen Physik, ed. by J. Nitsch et 
al, Bibliograph. Inst. Mannheim, Wien, Ziirich (1981) 65-84 (in German). 

[16] J. Ehlers, T. Buchert, Gen. Rel. Grav. 29 (1997) 733. 

[17] A. Friedmann, Generalized functions and partial differential equations. 
Prentice Hall, Englwood Cliffs (1963). 



14 



O. Heckmann, E. Schiicking, Zeitschrift fiir Astrophysik 38 (1955) 95 (in 
German). 

O. Heckmann, E. Schiicking, Zeitschrift fiir Astrophysik 40 (1956) 81 (in 
German). 

O. Heckmann, E. Schiicking, Encyclopedia of Physics 53 (1959) 489, Springer 
(in German). 

E.A. Novikov, Sov. Phys. JETP 30 (1970) 512. 

P.J.E. Peebles, The Large-scale Structure of the Universe, Princeton Univ. 
Press (1980). 

V. Sahni, P. Coles, Phys. Rep. 262 (1995) 1. 

J. Scrrin, Encyclopedia of Physics VIII. 1 (1959), Springer. 

A. Silbergleit, J. Math. Phys. 36 (1995) 847. 

Ya. B. Zel'dovich, Astron. Astrophys. 5 (1970) 84. 



15 



